Evidence for universal scaling in the spin-glass phase 
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We perform Monte Carlo simulations of Ising spin-glass models in three and four dimensions, as 
well as of Migdal-Kadanoff spin glasses on a hierarchical lattice. Our results show strong evidence 
for universal scaling in the spin-glass phase in all three models. Not only does this allow for a clean 
way to compare results obtained from different coupling distributions, it also suggests that a so far 
elusive renormalization group approach within the spin-glass phase may actually be feasible. 
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The characterization of the spin-glass phase in finite 
space dimensions remains as one of the prominent unre- 
solved problems in the physics of disordered systems. Al- 
though the spectrum of theories that have beenproposed 
to describe spin glasses is broad [H, 0, 0, 0, @, U 0,11, H , 
the main theoretical pictures are the simple scaling ap- 
proach (droplet picture) 0, S| and the replica symmetry 
breaking (RSB) scenario inspired by mean-field theory 
Because frustration and disorder are inherent ingre- 
dients of spin glasses, progress via field-theoretical calcu- 
lations has been difficult, at least below six dimensions 
. Thus most of the progress in the field relies on nu- 
merical studies that are also strongly limited mainly due 
to the numerical complexity of spin glasses. In fact, the 
difficulties are such that there is no good numerical evi- 
dence of whether a renormalization group (RG) approach 
may work in the spin-glass phase. 

We address some of the necessary features that one 
should observe in the spin-glass phase such that scaling 
theory and an RG approach may be potentially success- 
ful. To check for the applicability of a scaling theory 
within the spin-glass phase we assume a priori that a scal- 
ing approach works and check a posteriori whether our 
results are consistent with the scaling assumptions made. 
For this purpose we study finite-size scaling functions 
where we compare the behavior of several observables on 
a change of scale: We rescale the system size and plot 
observables as a function of a phenomenolo gica l coupling 
which in our case is the Binder cumulant We find 

that this procedure leads to scaling functions consistent 
with a universal scaling behavior in the spin-glass phase 
for all models studied: the Migdal-Kadanoff (MK) spin 
glass on a (three-dimensional) hierarchical lattice [12| , as 
well as the three- (3D) and the four-dimensional (4D) 
Edwards- Anderson (EA) short-range Ising spin glasses 
[3] ■ Of particular interest is the scaling function of the 
Binder cumulant itself since it provides a compact way 
to look at the complete RG flow. 

The MK spin glass serves as an example to illustrate 
the behavior of the different finite-size scaling functions 
in a simple scaling theory. This, in turn, follows from the 
possibility to exactly solve the model (numerically) using 



an RG decimation transformation. 

The physically interesting case of the 3D EA model 
is also the most difficult as the lower critical dimension 
is close [H|, i.e., that the spin-glass phase in this case 
is rather marginal. Hence we also study the model in 
4D. We observe similarities as well as clear differences 
between the MK and EA models. In particular, for all 
studied models we find clear evidence for the emergence 
of a universal scaling behavior in the spin-glass phase in 
the thermodynamic limit. Although our findings might 
be taken as a hint for the correctness of a simple scaling 
approach for the EA model, it is fair to observe that 
for very low temperatures we find an effective stiffness 
exponent 9 which is compatible with zero, as expected in 
the RSB or TNT @, H scenarios. Furthermore, for the 
system sizes studied, the fractal dimension of the surface 
of the low-energy excitations seems not to be equal to 
the space dimension d, in agreement with the TNT or 
droplet scenarios. We therefore lay the foundation for a 
simple approach that should allow future studies to check 
whether the disagreement with the traditional pictures 
is due to scaling corrections, or whether new theoretical 
descriptions are needed. 

Finite-size scaling approach. — For any observable 
0(L,T) as a function of the temperature T and the sys- 
tem size L, and the finite-size correlation length £(L,T) 
finite-size scaling (FSS) theory predicts [lj| that 

0(sL,T)/0{L,T) = F [f(L,T)/L;s] + (corrections), 

where s is a scaling factor. The corrections vanish in 
the thermodynamic limit. Because even the correlation 
length divided by the system size diverges in the spin- 
glass phase, this definition is inconvenient to examine the 
behavior of FSS functions within the spin-glass phase. 
Hence we use an alternative phenomenological coupling 
as the scaling variable which allows for a better visualiza- 
tion of the scaling functions in the spin-glass phase. We 
find that the Binder cumulant g(L, T) [see Eq. ([2]) below] 
works best since it is bounded in the interval [0, 1], thus 
leading to a compact picture of the scaling properties in 
the whole spin-glass phase. In the following we study 
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different FSS functions 

0{sL,T)/0{L,T) = F [g{L,T);s] + (corrections) (1) 

denned as a function of g(L,T). Next, we verify numer- 
ically whether it is sensible to study such FSS functions 
in the spin-glass phase and discuss our findings. 

Models and observables. — We consider the spin-glass 
Hamiltonian H. = — Yljj JijO~iO~j, where the sum is over 
all nearest neighbor pairs on the lattice and Cj £ {±1} 
are Ising spins. The couplings Jy are drawn from either a 
Gaussian (G), bimodal (B), link-diluted bimodal (D), or 
irrational (I) distribution which corresponds to a bimodal 
distribution where half of the (randomly chosen) bonds 
are multiplied by an irrational constant ci = (1 + v / 5)/2. 
The Hamiltonian is studied both on a hierarchical lattice 
with an effective space dimension of three and on a simple 
hypcrcubic lattice in three and four space dimensions. 
We use periodic boundary conditions for the EA model 
on the hypercubic lattices and free boundary conditions 
for the hierarchical lattice [l6|. The space dimension of 
the hierarchical lattice with G generations is d = 1 + 
ln(6)/ln(s), where b is the number of parallel branches 
and s is the number of bonds in series (we set s = 2 and 
b = 4 to obtain an effective space dimension of 3). The 
size of the system is L = s G . 

The spin-glass order parameter is given by q = 
(1/N) J^i °~i T i-, where a and r are two replicas of the sys- 
tem with the same disorder. For the hierarchical lattice 
we use [13] q = l/(2N L ) ( a * T i + a 3 T j) where the 
sum runs over all links Nl of the lattice. The spin-glass 
susceptibility \ is defined via x(L, T) = N[(q 2 )] av , where 
( • ) represents a thermal average and [ • ] av a disorder av- 
erage. The Binder cumulant g is defined as: 



g(L,T) 
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Finally, we also 



study the link-overlap qi = 
j, where the sum is over all links. 
Within the TNT picture the fractal dimension d s of 
large-scale excitations can be extracted from the variance 
of the link overlap a 2 ql (L,T) = [(qf) - ( qi )% v ~ 
where /i; = 9 + 2(d — d s ) varies with temperature [l8j |. 
The stiffness exponent 9 follows from the temperature 
derivative of the finite-size scaling function of the Binder 
cumulant Fq t9 via the quotient method [l9| : 



l + g*d g F g (g,s)\ 



g=g* 



corrections, 



(3) 



where g* is the value of the Binder cumulant at any given 
strong-coupling fixed-point [2(| and s is the scale factor 
used in the definition of the scaling function. Note that 
cumulants of the order parameter as well as the correla- 
tion length £/L are RG invariant quantities often referred 
to as "phcnomcnological couplings." In contrast to tra- 
ditional scaling analyses of the spin-glass phase which 



TABLE I: Parameters for the simulations of the 3D model 
with Gaussian (3DG) and irrational (3DI) disorder, as well as 
the 4D model with Gaussian (4DG) and bond-diluted (4DD) 
disorder. L is the system size, iV S a is the number of disorder 
realizations, N sw is the number of equilibration and measure- 
ment sweeps, T m i n is the lowest temperature and iV r the num- 
ber temperatures used in the exchange Monte Carlo method. 
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use bare (unrenormalized) couplings (e.g., temperature), 
here we use renormalized couplings (e.g., Binder ratio) 
as scaling variables, thus presenting a cleaner analysis. 

Computational details. — For the hierarchical lattice 
we use a similar procedure to Ref. 



2l|. An alternative 



way to calculate the link- and spin-overlaps on the hier- 
archical lattice is given in Ref. [22|. For L < 16 (G < 4) 
we use 7V sa = 5 • 10 5 samples, for L = 32 (G = 5) at least 
10 5 samples, for L = 64 (G = 6) at least 9 • 10 4 samples 
and for L = 128 (G = 7) at least 1.5 • 10 4 samples. All 
data are averaged over 10 3 independent configurations. 

For the regular lattices the simulations are performed 



using exchange Monte Carlo [23j, |24|. For systems 
with Gaussian disorder we use the equilibration test 
of Ref. [3], whereas for the link-diluted case we per- 
form a logarithmic binning of the data. Once the last 
three bins agree within error bars the system is equi- 
librated. For the 3D model with irrational [Gaussian] 
couplings T c = 1.47(3) [T c = 0.951(9)]. For the 4D link- 
diluted model with 35% dilution [25j [Gaussian disorder] 
T c = 1.0385(25) [T c = 1.805(10)]. For details see Table 
III 
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Finite-size scaling functions. — Wc first address the 
scaling function of the Binder cumulant F g for all mod- 
els studied and then discuss further scaling functions. In 
Fig. [T] we show our results for the FSS function F g for 
different models. F g displays in a compact way the RG 
flow of the Binder cumulant. Panels (a) and (b) show 
results for the hierarchical model where the scaling ap- 
proach is known to work. Panel (a) shows a comparison 
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FIG. 1: (Color online) Comparison of the finite-size scaling 
function of the Binder cumulant for different models, (a) 
Gaussian (open/black symbols) and irrational (full/red sym- 
bols) disorder for the hierarchical lattice, (b) Hierarchical 
lattice with bimodal disorder. The data converge slowly to 
the Gaussian limiting case (solid curve), (c) 3D EA model 
with Gaussian (open/black symbols) and irrational (full/red 
symbols) disorder, (d) 4D EA model with both link-diluted 
(full/red symbols) and Gaussian (open/black symbols) disor- 
der. The dashed lines in all panels represent a droplet scaling 
behavior with 9 = 0.27 for the hierarchical, 9 = 0.2 for the 
3D and 9 = 0.75 for the 4D models, respectively. 



between Gaussian and irrational coupling distributions, 
whereas panel (b) shows how the FSS scaling function for 
the bimodal coupling distribution — amidst strong finite 
size effects — converges (slowly) towards the FSS function 
of the Gaussian (and irrational) cases. The convergence 
is slowest close to g = 1 (T = 0) because for low T en- 
tropic effects become relevant [26[. These results are a 
clear indication for universal FSS in the spin-glass phase 
for the MK model. Panel (c) shows a comparison of the 
Gaussian and the irrational coupling distribution in the 
3D EA model. Finally, panel (d) shows a comparison 
between the Gaussian and the link-diluted bimodal EA 
model in 4D. The results for the 3D/4D EA models are 
consistent with the scaling hypothesis and indicate a uni- 
versal scaling behavior. The broken lines indicates how 
the FSS function should depart from g = 1 assuming 
that simple droplet scaling is correct, i.e., that the slope 
is given by the exponent 9 through Eq. ([3]). While this 
works perfectly for the MK model, there is a clear differ- 
ence in the case of the EA model. 

In order to show that the scaling collapse in Fig. Q] is 
not coincidental we show in Fig. [2] the FSS function of 
the spin-glass susceptibility. The top panel shows a com- 
parison of Gaussian and irrational disorder for the MK 
model, whereas the bottom panel shows a comparison for 
the 3D EA model. Again, the data show strong evidence 
of a universal scaling behavior. The data for the MK 
model seem to fall onto a straight line which underlines 
the very simple scaling behavior in this model. 

In Fig. [3] we show results for the link-overlap variance. 
The top panel shows data for the MK model with Gaus- 
sian and irrational disorder distributions. In contrast to 
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FIG. 2: (Color online) Comparison of the finite-size scaling 
function of the susceptibility \ f° r the hierarchical lattice 
(top) with the 3D EA model (bottom). In both models the 
results for the Gaussian (open/black symbols) and the irra- 
tional (full/red symbols) coupling distribution are shown. 
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the other FSS functions discussed so far we find sizable 
scaling corrections. The data are less conclusive but there 
is a clear trend that the curves for the two different dis- 
tributions converge to a single master curve. The bottom 
panel shows data for the aforementioned disorder distri- 
butions in the 3D EA model. The data are similar to the 
MK case and show small corrections to scaling. The bro- 
ken lines (or the region between them) correspond(s) to 
the expected value of the scaling function at g = 1 in the 
thermodynamic limit under droplet scaling assumptions 
(here we use 6 = 0.20(5) and d - d s = 0.42(3) Q in the 
3D EA case; 9 = 0.27 and d s = d - 1 = 2 in the MK 
case). Neither the droplet nor the RSB picture (where 
F a 2 — > 1) extrapolate to a consistent value for F a i in 
the limit g — > g(L, T = 0), in agreement with Refs. 

Summary and discussion. — Studying the behavior 
of several FSS functions we have found evidence for a 
universal scaling behavior in the spin-glass phase. The 
existence of these universal FSS functions allows us to 
perform for the first time a precise comparison of re- 
sults in the spin-glass phase between different coupling 
distributions. We find that neither the simple droplet 
nor the RSB picture extrapolate to consistent values for 
the scaling functions in the zero-temperature limit. For 
both pictures it can be argued that these inconsistencies 
might be due to scaling corrections 2?], 28[. Our results 
suggest that in such a case not only the scaling behav- 
ior in the thermodynamic limit is universal, but also the 
leading scaling corrections. Although our results do not 
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FIG. 

function of the link-overlap variance a% for the hierarchical 



3: (Color online) Comparison of the finite-size scaling 
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model (top) and the 3D EA model (bottom). For both MK an 
EA models the results for the Gaussian (open/black symbols) 
and the irrational disorder (full/red symbols) are shown. The 
broken line (the range between the broken lines, respectively) 
correspond(s) to the expected value of the scaling function for 
g — * 1 in the thermodynamic limit within the droplet model. 



allow us to discriminate between the different scenarios 
that have been proposed for the nature of the spin-glass 
phase, they allow for a new, parameter-free way of look- 
ing at the problem. 
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